Effects of Eye-phase in DNA unzipping 
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Abstract 

The onset of an "eye-phase" and its role during the DNA unzipping is studied when a force is apphed 
to the interior of the chain. The directionahty of the hydrogen bond introduced here shows oscillations in 
force-extension curve similar to a "saw-tooth" kind of oscillations seen in the protein unfolding experiments. 
The effects of intermediates (hairpins) and stacking energies on the melting profile have also been discussed. 
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Molecular interactions play a key role in living organisms. Recent advances in experimental 
techniques have allowed nanomanipulation in single biological molecule and made possible to mea- 
sure these interactions [ij. The aim is to exert a force in the pN range by optical tweezers, atomic 
force microscopyjetc and characterize the molecular, elastic, structural and functional properties 



of bio- molecules [2, l3|] • In typical experiments of double stranded DNA (dsDNA) unzipping, a force 
is applied to the ends of the chain (Fig. la) and one studies the force-extension curve which shows 
the elastic properties and the gross features that reflect the local 'G-C vs 'A-T' content along 



the sequence 



mm. 



Moreover, the force-temperature diagram shows, below the melting 



temperature, that the over-stretching force [s] decreases nearly linearly with temperature. 

Theoretically, DNA unzipping may be studied in different ensembles depending on the ex- 
periments. For example, Atomic force microscopes (AFM) work in constant distance ensemble 
(CDE) while magnetic bead uses the concept of constant force ensemble (CFE). The prediction of 
unzipping transition based on interacting Gaussian chains [10] raised a lot interest and now results 



are available from dynamical approach [111] , exact solutions of lattice models 



12 



simple models 

n 



of quenched-averaged DNA [ij, [l5| , numerical simulations and scaling analysis [16|, llTI] . Recently 
for a model of interacting polymers where any monomer of one chain can interact with any monomer 
of the other chain (we call it model A), the role of an intermediate entropy-stabilized phase was 



recognized and a force-induced triple point 



18 



In most of the models studied for dsDNA 



in a force-tem per ature plane was established. 

Q, Q, Q, QQ, 



12 



a monomer i of one strand 



can only interact with the i-th monomer of the other strand, which is similar to the models of 

□ 

DNA (we call it model B) proposed earlier by Poland and Scheraga [20|]. These models do not 
take into account the directional nature of the hydrogen bond and underestimate the entropy by 
restricting the formation of hydrogen bonds. Thus these models may only give a limited picture of 
the unzipping transition and do not allow to study the effect of intermediate states 21 ] . 



So far, most of the theories of DNA unzipping have focused on the thermodynamic limit, and 
therefore, consist of a few parameters typically in the form of an effective base-pair interaction 
in a simplified way. In all the single molecule experiments, a finite size chain is used and hence 
no "true phase transition" can in principle be observed. Still, the "phase transition" observed in 
such experiments may be considered as real if the length of the chain exceeds the characteristic 
correlation lengths. It is now becoming possible to go to the other limit of studying shorter 



22 



231 ]. The purpose of this letter is to provide 



segments (~ 10 base pair) at coarse grained level 
exact results of a semi-microscopic model of short chains by incorporating the directional nature 
of hydrogen bonds and then propose a method to study the effect of molecular interactions right 
at the individual base pair level and their role on melting profile. 

In the following, we adopt a more realistic model of DNA, which may be defined in any dimension 




FIG. 1: The schematic representation of DNA unzipping by force (g) apphed at (a) one end of the two 
strands; (b) at the interior of the strands. Figs, (c-e) represent the possible conformations of model C. Here 
(c) represents the completely zipped state, (d) a non pairing configuration and (e) a partial bound state 
with hairpin loop. 



2j] (herein after we call it model C). A similar model has been used in Ref. 24 in the context 
of relative stabilities of DNA hairpin structures. The model takes care of important shortcomings 
of model B and also incorporates some additional features like existence of intermediate states, 
effects of stacking energy, excluded volume properties of nucleotides and the directional nature of 
hydrogen bonds. 

Although the importance of bubble formation during thermal melting has been recognized, no 
attempt has been made so far to experimentally explore the phase diagram when a force is applied 
to the interior of the chain as shown in Fig lb. Such situations occur in many biological processes, 
for example, during gene-expression, RNA forms bubbles or "eye-type" conformations on DNA. 
Therefore, we consider two cases: (i) force has been applied at the end of the chain (Fig. la, 'END' 
case or Y case), and (ii) at the middle of the chain (Fig. lb, 'MID' case). The contribution to 
energy by this force, g, is —2gx, where 2x is the absolute distance in the x-direction between the 
two strands at the point of application of the force. 

We model the two strands {e.g. A-T) of a homopolymer DNA by two mutually- attracting- 
self-avoiding walks (MASAWs) on a square lattice as shown in Figs, l(c-e). The bases are 
associated with the link between two monomers of a chain as depicted in Fig. 1. In one strand the 
bases point towards the right while on the other they are on the left, as one traverses the chains 
sequentially. We associate a contact energy — e (effective base pair interaction) with each pairing 
between complementary strands only if the bases are nearest neighbors (short range nature of the 
hydrogen bond) and approach each other directly without the strands coming in between [Fig. Ic]. 
Fig. le shows the possibility of formation of hairpin (which is not possible in model B) in a single 
strand of DNA. However, in this case, non-native contribution has been taken into account but no 
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FIG. 2: Variation of critical force {gc{T)) (END and MID case) as a function of temperature (T) for model 
B and C. 

apparent weight has been assigned to stem as it is made up of same nucleotides. 

The partitions function (^at) of the system under consideration can be written as a sum over 
all possible configurations of MASAWs i.e. J2m x C{m, x) exp(/3me) exp{f3gx), where /? = l/ksT 
is the inverse temperature, ks being the Boltzmann constant. C{m,x) is the number of distinct 
conformations walks of length 2N having m number of intact base pairs whose end (or mid) points 
are at a distance x apart. We have obtained C{m, x) for < 16 and analyzed the partition function 



using exact enumeration and series analysis technique 



25 



26l |. We prefer this technique because 



it can predict various phases of the system [27] quite effectively and the scaling corrections can 



251, 



26l |. To achieve the same 



be correctly taken into account by a suitable extrapolation methods 
accuracy in Monte Carlo, a chain of two orders of magnitude larger than in the exact enumeration 
method is required 



28l |. We set e/ks = 1 and calculate the reduced free energy per base pair from 
the relation G(T,g) = lirn^^oo ji^og Z{T, g) = log fj,(T,g) 26(]. The limit ^ oo is achieved by 



using the ratio method [25| for extrapolation. The transition point can be obtained from the plot 
of G{T, g) versus T or from the peak value of g^fy^p-- 

A force-temperature [g — T) phase diagram of model B and C for end and mid case is shown 
in Fig. 2. The qualitative features of the phase diagram obtained here may be compared with 
experiments [§]. The phase boundary separates the zipped and the unzipped state. At T = 0, the 
critical force can be found from a simple analytic calculation and is equal to 0.5 [isl that is in 
agreement with the one from Fig. 2. 

Experimentally, the melting profiles are obtained by monitoring the change in the UV ab- 
sorbence with temperature which provides the information about the fraction of open base pairs 
and the melting temperature is defined when half of the total base pairs get open [29|]. Another 



quantity of experimental interest is to monitor and measure the end separation by varying the force. 
We calculate these quantities from the expressions < ni >= J2nT'C{m,x)eyip{l3me)exp((3gx)/ZN 
and < X >= ^xC{m,x)exp(l3me)exp{Pgx)/Zi\i and plot their variations with temperature and 
force in Figs. 3 and 4, respectively. It is evident from Fig. 3 that, with increasing temperature 
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FIG. 3: Variation of scaled < to > (a) with temperature (T) at constant force g = 0.25 (END) and 0.50 
(MID); (b) with force {g) at constant T = 0.3. Here sohd and dashed hnes represent the END case, while 
dot and dot-dashed line represent the MID case for model B and C respectively. 
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FIG. 4: Plot of scaled < a; > (a) with force {g) at constant T — 0.3 and (b) with temperature (T) at 
constant g — 0.25 (END) and 0.5 (MID). The lines have same meaning as of Fig. 3. 

(or force), the number of intact base pairs decreases and there is a sharp transition from a zipped 
state to the unzipped state. 

Remarkably, to break the same amount of base pairs at low temperature one requires almost 
double the force if it is applied in the middle of the chain, rather than the end (Fig. 3b) consistent 



with the exact results 13|]. Near the melting temperature, fluctuation dominates and less than 
double the force is required for unzipping the chain from the middle. 

Another interesting observation is the variation of average elongation with force which shows a 
monotonically increasing trend (Fig. 4a) at constant temperature and approaches unity. Meanwhile 
the variation of extension at constant force shows a sharp rise with temperature (Fig. 4b) and then 
a slight decline to approach a value below unity. At constant temperature, there is a transition 
from the zipped to the unzipped ("rod-like") state. Keeping force constant, when temperature 
is varied, there is a transition and the chain acquires conformations close to the rod-like states. 
As this temperature is still low, with further increase in temperature, the entropy of the system 
increases and the chain acquires coil-like state at higher temperature and thus average distance 
decreases. 

We also study the probability distribution curves P{x) with x and P{m) with m for 
model B and C defined by P{x) = YyTn^NirriTx) exp{l3me) exp{Pgx)/Z]\f and P(rri) = 
J2xC^NiiT^, x) exp{f3me) ex.p{Pgx)/ZN respectively. In Figs. 5(a-f) we have shown P{x) and P{m) 
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FIG. 5: Figs, (a-f) represent the P(x) and P{m) of END case (a-c) and MID case (d-f) for different forces 
at constant T=:0.3. Fig. 5e shows the signature of "eye-phase" of even widths in the form of oscillation for 
the mid case in model C but absent in model B. 

for different values of force and a fixed temperature T = 0.3. The x-component of the distribution 
function gives information about the states of dsDNA. The maxima of P{x) at x correspond 
to the zipped state for a given set with g = 0.1 and T = 0.3. For small forces, thermal fluctuations 
are too weak to unzip the strand and the DNA remains in the zipped state. This is being reflected 
in the structure-less distribution function with a well defined peak at the most likely value of the 
extension (Figs. 5a and 5d). For both the models at higher force and at the same temperature 
{g = 0.6 for end and 0.95 for mid respectively at T = 0.3), thermal fluctuations have been sup- 
pressed by a strong force and strands are found in the segregated "rod-like" state with more or less 
an identical distribution peaked at maximum extension. However, striking differences are observed 
in the probability distribution curves for model B and C at intermediate forces below transition 
line {g — T plane) when the force is applied in the middle. No such differences are observed for 
the end case. For the mid case, the probability distribution curve of model C shows strong oscil- 
lations, whereas for model B no such oscillations are observed. This indicates that the model B 
may be described by two state models. However, model C shows that certain intermediates states 
("eye-shape" of even width) are more favorable than others. 

It is interesting to note that the g — T phase diagram of model B and C are almost identical 
without any re-entrance at low temperature. This may be because of the fact that the energy of 
the ground state and the unbound state for both models are the same. At the center point of 
the 'Y' (end case), the phase boundary is determined by a balance of the net force -2gx and the 
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FIG. 6: Schematic diagrams showing the different conformations keeping 2x (< iV) distance constant. For 
model B, (b) and (c) has same weight while for model C, (c) has more weight than (b). 
unzipping potential —em with associated entropy. In the thermodynamic limit, directionality of 
the hydrogen bonds, which appear in the form of entropy of the partial bound states does not play 
a crucial role in this balance. The absence of re-entrance is understood with the zero entropy of 
the ground state for both models. As shown recently by Kapri et al. [13] for mid case, in the 
CDE there is a possibility of a coexistence region that is better thought of as an "eye-phase" = 
two 'Y' joined together. In this case, the separation at the point of application of force is smaller 
than the fully open case and thus such conformations statistically have more weight than the other 
conformations. 

Since the bottom end is kept fixed, the top side of the strand may open due to thermal fluc- 
tuations (Fig. 6b) and form a partial bound state as shown in Fig. 6c which is more stable than 
Fig. 6b. Therefore, in model C, the half of the chain undergoes an unzipping transition while 
other half due to the non-native contacts, shows the combined effects of unzipping and slippage 
(shearing) transition. Thus in the model C, the transition appears more smoother than the model 
B which can be seen in Fig. 3b. The consecutive peaks in P{x) vs x curve (Fig. 5e) represent 
the "eye-phase" of even widths which contribute most to the partition function. Because of native 
contacts, the contribution of the "eye-shape" conformation is significantly less and small thermal 
fluctuations are sufficient to unzip the chain in model B. 

We substantiate our arguments by extending calculation in CDE also. The partition function 
in CDE may be defined as Z_m{x, T) = Y^m exp(/3me). The two ensemble are related by Zn{T^ g) = 
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14( 1 . The free energy is given by the relation F]\f {x,T) = — T In {x,T). 



In CFE the average separation < a; > fluctuates while in CDE one measures the average force 
to keep the separation constant given by the expression < g >= ^^^g^''^^ at constant temperature 
^j. The force-extension curve thus obtained is shown in Fig. 7 for model C which also shows 
oscillations for the 'MID' case but constant for the 'END' case. Though a somewhat similar effect 
has been seen experimentally in molecules like Titin |3Q], DNA has not been probed so far. The 
average of force obtained here (0.2 ±0.05 for 'END' case and 0.35 ±0.05 for 'MID' case at T = 0.5) 
also matches with the value shown in Fig. 2. 

The stacking energy in case of homopolymers gets adsorbed in the effective base pair interaction. 
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FIG. 7: Plot of < ,g > with distance (x) at constant T = 0.5. 
To see this, we associate an additional energy between two consecutive parallel base pairs (only 
possible in model C) and found that there is no qualitative change in the phase diagram except 
shift in the transition temperature. If stacking energy is negative (attractive interaction) chain 
becomes stiffer. However, it does not change the nature of oscillations as observed in probability 
distribution curves for the 'MID' case. 

The exact results on short chains of a new semi-microscopic dsDNA that incorporates the 
directional nature of hydrogen bond show unequivocal signature of an "eye phase" , without going 
to the long chain limit. This happens for the case with a force acting in the middle of the dsDNA. 
The variation of elongation due to the force in different ensembles has different behaviour both 
qualitatively and quantitatively. We anticipate that refinements in high precision single molecule 
experiments will be able to verify these predictions. 
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